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We perform the first bottom up analysis of an idealized ecosystem using Buckingham's Pi theorem. 
Our analysis requires only three assumptions, (i) that the ecosystem is in a dynamically balanced 
steady state, (ii) that some (non- unique) method exists that can differentiate different functional 
groups of individuals within the ecosystem and (iii) that these functional groups are connected to 
each other by the flow of resource. The expression that we obtain relates the level of complexity that 
the ecosystem can support to intrinsic macroscopic variables such as density, diversity and charac- 
teristic length scales for foraging or dispersal, and extrinsic macroscopic variables such as habitat 
size and the rate of supply of resource. This expression relates these macroscopic variables to each 
other, generating commonly observed macroecological patterns. Since we do not require any details 
of the ecosystem function to obtain these broad macroecological patterns this may explain why 
they are ubiquitous. Departures from our relationship may indicate that the ecosystem is in a state 
of rapid change, i.e., abundance or diversity explosion or collapse. Our result provides normalised 
variables that can be used to isolate the trend in one ecosystem variable from another, providing a 
new method for isolating macroecological patterns in data. Quite generally, a dimensionless control 
parameter for ecosystem complexity emerges from our dimensional analysis. This, we suggest, will 
be a control parameter in dynamical models for ecosystems based on energy flow and conservation 
and will order the emergent behaviour of these models. 

PACS numbers: 87.23.-n,89.75.Fb 



I. INTRODUCTION 

The rapid increase in the availability of large-scale 
ecological data, and the recent rise of the discipline of 
macroecology jll, Hj have simultaneously increased our 
knowledge of global patterns and stimulated the search 
for the underlying processes that determine them (see 
3). Examples of the large-scale, macroecological, pat- 
terns Q to have emerged from empirical analyses include 
the species-area and the species-latitude relationships, 
and trends in density and diversity with body size. These 
are broad scale patterns and generalised rules rather than 
mechanistic processes, and a range of theories have been 
proposed to explain why these patterns emerge. This has 
led to the development of a perspective in which the de- 
tailed biological characteristics of the species (traits) do 
not determine their abundance and that the processes 
affectirig community structure can be considered as neu- 
tral 043, see also [8, 9]. The fact that these patterns are 
both approximate and ubiquitous, suggests that they do 
not reflect the details of how ecosystems operate, rather 
that they emerge from underlying general constraints. 

Rather than exploiting statistical constraints, as in 
maximum entropy inference [gI. [lol - [l^ . we focus on physi- 
cal and dimensional constraints and how these relate the 
level of complexity that the ecosystem can support, to 
intrinsic variables such as density, diversity and charac- 
teristic lengthscales for foraging or dispersal, and extrin- 



* Electronic address: 's.c.chapman@warwick.ac.uk| 



sic variables such as habitat size and the rate of supply 
of resource. Any model or function that relates physical 
quantities has a fundamental property; that the interde- 
pendencies which it captures must be independent of the 
units in which the quantities are measured, that is, it 
must be exp ressible completely in dimensionless parame- 
ters |13l-[l5l|. As a consequence, if we observe a particular 
set of ecosystem variables, since these variables possess 
physical dimensions (mass, length, time) there are di- 
mensional constrains on the relationships that can exist 
between them. Importantly, formal dimensional analy- 
sis yields these dimensional constraints without detailed 
knowledge of the functioning of the ecosystem; it does 
not specify a unique relationship between the variables. 
We can obtain this if in addition there are known physical 
constraints, such as conservation principles. The physi- 
cal constraint that will be central to this paper is that 
of an ecosystem which adapts dynamically to changes in 
external parameters to maintain a balance between the 
rate of uptake and of utilization of resources taken over 
the ecosystem as a whole. The concept of ecosystems 
which compensate dynamically to remain in homeostasis 
has been explored by [H, [ijj but not developed using 
dimensional analysis. 

In this paper we will only require three key assump- 
tions, (i) that the ecosystem is in a dynamically balanced 
steady state, (ii) that some (non- unique) method of ob- 
servation can differentiate different functional groups of 
individuals within the ecosystem and (iii) that these func- 
tional groups are connected to each other by the flow 
of resource. This will yield an expression that relates 
the level of complexity in the ecosystem as quantified 
by (ii) to its characteristic intrinsic and extrinsic vari- 
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ables. This expression constrains the relationships that 
can exist between these (dimensional) ecosystem vari- 
ables which we will see is reflected in overall macroe- 
cological trends that are observed. Further, we obtain 
the relationship between these trends. At minimum this 
leads to the dimensionless, or normalised variables that 
need to be constructed to isolate variation in one ecosys- 
tem variable from another; this provides a new method 
for isolating macroecological trends. At maximum, we 
obtain a dimensionless control parameter and relate it 
to the level of complexity that the ecosystem can sup- 
port, which points towards 'thresholds for life', that is, a 
parametrization of the minimum level of complexity that 
a dynamically balanced ecosystem can support. 

Process based models of ecosystems that aim to predict 
macroecological patterns tend to fall into two approaches. 
The first of these relates area, diversity and abundance 
by means of models for dispersal, occupancy and coexis- 
tence /conipetition and speciation in physical space (see 
e.g. [3, S E, [isl)- The second relates abundance, body 
mass and metabolic rate by means of models and con- 
straints for the availability and fiow of resource ([l^l, see 
the review of [21|] ) . The relationship that we obtain here 
links the key parameters of both these approaches, sug- 
gesting a synthesis of them. 



II. FORMAL DIMENSIONAL ANALYSIS AND A 
BOTTOM UP- APPROACH TO AN ECOSYSTEM 

We will build our understanding by first considering 
the simplest possible idealized ecosystem and then suc- 
cessively increasing the level of complexity; at each stage, 
additional variables are introduced. Our approach is to 
use formal dimensional analysis at each stage to find 
the constraints that act on this general description of an 
ecosystem. Importantly, we seek to describe an 'observed 
macroscopic ecosystem', that is, the variables that we will 
ultimately identify include observed intrinsic properties 
such as density, diversity, body size, and metabolic rate. 
Introducing progressively more specific detail inevitably 
introduces more characteristic variables- this procedure 
could be taken further to explore specific detailed ecosys- 
tem models by the input of more detailed ecosystem func- 
tional properties into the dimensional analysis. Our aim 
here is rather to explore the ecosystem constraints that 
emerge for the minimum set of assumptions and model 
inputs. 

Any model for a physical system that depends upon V 
variables Qi,Q2t-Qv must provide a function F that 
relates them: F(Qi, Q27 --Qv) = C where C is some 
dimensionless constant [15|. The physical system that 
we discuss here is that which captures general physi- 
cal aspects of ecosystem function, specifically the uptake 
and utilization of resource. The essential idea of for- 
mal dimensional analysis is that this function can only 
depend upon dimensionless parameters ni,n2,..nM so 
that F = F(ni,n2, ..IIm) only. These dimensionless pa- 



rameters TIi..m{Qi..v) are formed directly from the vari- 
ables Qi,Q2, --Qv- Formal dimensional analysis as in 
Buckingham's 11 theorem ( [T3l - [l5| , see also [13) is simply 
the process to obtain these dimensionless groups of vari- 
ables. If one then has additional information about the 
system, such as a conservation property, the Tli..MiQi..v) 
can be related to each other to make F explicit. If the 
V relevant variables Qi..v are expressed in W physical 
dimensions (i.e. mass, length, time) then from [13| there 
will be at least M — V — W dimensionless parameters or 
groups ni..M(Qi..y) which we now progressively identify. 

A. One kind of uniformly distributed single cell 
organism. 

The ecosystem is composed entirely of single cell organ- 
isms that are of uniform type: they have the same func- 
tion and structure and same typical metabolic rate R, di- 
mensions [M][L]2[T]-3. They are uniformly distributed 
in a habitat of size L in D dimensions with density n, di- 
mensions [L]-^. The available resource (sunlight) is de- 
livered at rate P per unit area, dimensions [M] [T] ~^ and 
the (dimensionless) fraction a taken up by the ecosystem 
is a constant. There are 5 variables (i?, L, P, a, n) and 3 
physical dimensions (mass [M], length [L], time [T]) so 
we have 2 parameters (dimensionless groups) which are 



These are related by the physical constraint that the sys- 
tem is in dynamical balance so that the rate at which en- 
ergy is taken up over the ecosystem is the rate at which 
it is utilized, so that: 

^ ^ uL- (2) 
R ^ ' 

which is Hi =112. This expression simply tracks the flow 
of energy into and through the ecosystem- it assumes that 
all other processes necessary for the ecosystem to func- 
tion, such as the recycling of resources such as Nitrogen, 
occur. Introducing a typical metabolic rate for the single 
cell organisms has fixed an energetic minimum threshold 
for life which is Hi — 1, that is, one cell in the habitat 
(one cell ecosystem). 

B. More than one kind of single cell organism. 

We next consider single cell organisms that can be dif- 
ferentiated from each other, cither by their function or 
their structure, or both. These would represent distinct 
types or species and different methods for categorizing 
and distinguishing individuals will yield different sets of 
species. What will follow will be independent of the pre- 
cise details of this differentiation method, we only need 
assume that such a differentiation is now possible. In 
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the ecosystem there are S types (species) and there is 
a species label k — 1..S, the density of the fc*'' species 
is n{k). We can always define an average density of the 
single cell organisms: 



1 ^ 



(3) 



fe=i 



so that the variable n in ([T]) is now replaced by n and 
112 = nL^ . The additional variable S is dimensionless; 
so that we now have 6 variables and 3 dimensions, and 
so 3 dimensionless parameters with Ha = S. These are 
again related by the physical constraint that the system 
is in dynamical balance: 



aPL^ 
R 



s 

E 

fe=i 



n{k)L 
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fiSL' 



(4) 



Expression ^ now encapsulates the idea of primary 
producers- one or more of the species is responsible for 
the uptake of resource with efficiency a. The other 
species 'feed off' this primary producer either by graz- 
ing, predation or uptake of waste. 



C. Multicellular organisms 

We now consider more complex organisms that are 
multicellular. All the organisms live in an ecosystem and 
are connected to the primary producers by the flow of re- 
source, either directly or indirectly by predation, or both. 
It is now possible to distinguish types of organism and we 
will label the different types or categories distinguished 
in this way with index p. Again, the results to follow 
will not depend upon the precise details of how individ- 
uals are assigned to any of the p categories, simply that 
such an assignment can be made. Different methods for 
categorizing the individuals in an ecosystem [23j will or- 
ganize individuals into groups or categories of different 
p, this categorization may focus on the functional role 
of individuals in the ecosystem such as niche or trophic 
level, or may focus on stage of development, or other 
factors. Organisms falling into a given p*^ category or 
group will be clustered around an average body size, on 
average they will be composed of B{p) cells (this is typi- 
cally observed Q), and will have average metabolic rate 
RB{p), the per-cell metabolic rate R now corresponds to 
an ecosystem average over these multicellular organisms. 
There is a non- trivial correspondence between average 
size B{p) and how complex an organism can be. Within 
each p there will be a number k — l..S{p) of differen- 
tiable species each with density n(k,p) all with average 
body size B{p) and with average density, for that p of 



l{p) =< n{k,p) >k= 



1 



Sip) 



S{P) 
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Our relevant variables of density and diversity now re- 
late to the observation of a given category that is 
embedded in the ecosystem; we observe n* = n(p*), 
S"* = S{p^). Individuals in the observed category also 
has a characteristic average size B^, — B(j)^) (number 
of cells so dimensionless) which introduces an additional 
variable giving a total of 7 variables and hence 4 dimen- 
sionless parameters: 



R 



U2 = n,L^,U3 = S,,Ui = B, (6) 



The physical constraint of a dynamically balanced 
ecosystem is now: 



S(P) 



^ ^ n(fc, p)B{p)L^ = J2 n{p)S{p)B{pW 



R 

p fc=i p 

(7) 

We can write ([7]) in terms of our observed p''* category: 
-^^^^n.g.i..L^E^g^iyjg, (8) 



R 



n{p^)S{p^)B{p^) 



We then have an expression of the form: 



aPL^ 
R 



(9) 



fc=i 



where ^(p*) is a dimensionless function; l/(\I'(p*)) is 
the fraction of the total rate of resource supplied to the 
ecosystem that is utilized by the observed (pi^) category. 
Importantly, all of the species and categories are con- 
nected into the same resource flow, so that fundamen- 
tally, ([9]) will constrain how observed density, diversity 
and body size are related to each other across the ecosys- 
tem, with consequences for macroecological patterns as 
we will discuss. 

These dimensionless parameters ([B]), and their rela- 
tionship ([9]) express the following fundamental proper- 
ties of the idealized ecosystem. There is a building block 
on which life is organized- the single cell which has a 
definable typical metabolic rate, R. There is then the 
physical property of resource flow: that the single cell 
metabolic rate, along with the rate of uptake of resource 
to the ecosystem aPL^ constrains the number of cells 
the ecosystem can support, which is Hi. The detailed bi- 
ological and ecological properties of the ecosystem then 
organise these Hi cells into a complex network of species 
and groups of species, observationally these are charac- 
terized into functional units which have an average body 
size, density and diversity. Hence, the observed ecologi- 
cal variables are in a macroscopic sense related to each 
other by the physical property of resource flow. 



D. Non uniform distribution of individuals in 
space. 

Generally, organisms will not be uniformly distributed 
in space so that the density depends on the length scale 
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r over which it is observed, so that n = n{r,k,p) and 
similarly, the efficiency of the primary producers, which 
depends on their density, is a = a{r). The lengthscale 
over which the density varies can either arise from how 
individuals subdivide and grow, forage, or other forms 
of influence they have on each other and on the environ- 
ment. This is important since in (l7][9|) the density refers 
to that measured over the habitat of the entire ecosystem 
on scale L and any observation will be on a more local 
scale r << L, which in turn relates to the lengthscale 
of over which the density varies. This will introduce a 
variable for the scale of observation of the p*'* category; 
r* (dimension [L]) finally giving 8 variables so 5 dimen- 
sionless parameters: 



Hi 



a PL^ r 
, n2 - , Ha = , n4 - , Hs - ^ 

Jri L 

(10) 

where a* = a{r^,) so that all variables refer to a con- 
sistent set of observations on lengthscale r*. The ob- 
served density can be generally expressed as n{r,k,p) = 
n{L,k,p)/g{k^p,r/L) where g is dimensionless; similarly 
a(r) = a{L)/ga{r/L). If the categorization p is based 
on function and structure then one can anticipate that 
an average of g over the S{p) species in the category is 
meaningful so that 

^ Sip) 

n{r,p) =< n{r,k,p) >k= -^i—- V n{r,k,p) 



S{p) 



n{L, k,p) 



Sip) {^^gik,p,r/L) 



n{L,p)/gip,r/L) (11) 



With notation g^ = g{p^,) expression © becomes, for the 
observed : 



^^^^ ^Y.^(L'P)SiP)Bip)L^ 
p 

= J2 n{r,PmP, r/L)S{p)B{p)L^{U) 



= n S B g { — )L" V "(P)g(P' r/L)S{p)B{p) 
* * * * L ^ n{p^)g{p^,r^/L)S{p^)B{p^) 

or writing the spatial variation in a single function 

G{r/L)=gir/L)/g^{r/L) 

= n,S,B,G{^)L°^{p,) (13) 

Spatial variation in density thus leads to spatial trends 
in diversity, or species- area rules; this has arisen quite 
generally as a consequence of the physical constraint 
and we will discuss this next. 

Essentially, ([T^ is: 



Hi = n2n3n4G(n5)*(p*) 



(14) 



The dimensionless functions ^ and G capture all the de- 
tails of the ecosystem function and reflect how complex 
the ecosystem is. The level of complexity that an ecosys- 
tem can support increases with the number of different 
possible configurations that are possible and this is con- 
strained by the number of cells the ecosystem can sup- 
port. For example, if there is only one distinguishable 
kind of organism in the ecosystem, there is only one p 
and S value and Vl/ = 1, we essentially recover ([2]) where 
the average metabolic rate of the organisms is . The 
threshold for one such organism to be supported by the 
ecosystem is ni/n4 = 1 or a^,PL^ /(RB^,) = 1. As the 
ecosystem becomes more complex, ^I^ > 1 and each p cat- 
egory utilizes a smaller share of the total resource sup- 
plied. The dimensionless function ^I^ thus operates as an 
order parameter of the ecosystem which reflects the level 
of complexity. The corresponding control parameter is 
Hi . A physical analogy to this is the relationship between 
the Reynolds number in turbulence and the number of ex- 
cited modes or degrees of freedom. The Reynolds number 
is the ratio of a rate of energy input on the largest, driv- 
ing scale to a rate of energy dissipation on the smallest 
scale, as is Hi here, and dimensional analysis, along with 
the assumption of steady state (no energy pile up) is suf- 
ficient (see eg to constrain the numlDcr of degrees of 
freedom to grow with increasing Reynolds number. 



III. CONSTRAINTS ON MACROECOLOGICAL 
PATTERNS WITHIN AND ACROSS 
ECOSYSTEMS 

Observations both within and across ecosystems con- 
sist of specifying a method for classifying individuals into 
particular groups or categories and then for each of the 
p*'' categories, observing the average density, diversity, 
bodysize, and metabolic rate. From the constraint 
we see that these variables are not independent, and (jl3p 
suggests relationships between them which we will now 
discuss. 

Let us consider that a scheme for classifying individu- 
als is consistently adopted, and observations of average 
density, diversity, bodysize, and metabolic rate of these 
categories are made. These observations simultaneously 
collect a range of values of n* , 5* , B^ for a given , L 
and P. We will first consider the case where 'similar' 
ecosystems are compared, or where a comparison is made 
within a single ecosystem, that is, the order parameter ^ 
is not varying. Subsets of the variables in equation 
will then show functional relationships provided that the 
remaining variables do not strongly vary, this has been 
found for example by [13, HH] who demonstrate a rela- 
tionship between species richness, area and a measure of 
productivity. 

If one variable dominates our expression and 4* is held 
constant predicts the following patterns or trends: 

• Diversity and Wright's Rule: cx aPL^, so that 
the number of species (diversity) increases with 
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total net productivity integrated over the habitat 
rather than productivity alone; this is Wright's rule 
[26j . Whilst Wright's rule is to some extent eco- 
logically trivial (a greater net energy input allows 
more individuals, see [23|) the interesting aspect 
of this result is that it predicts an increase in di- 
versity (richness) and not just individuals. As we 
would anticipate from the relationship between 
productivity and species diversity also varies with 
spatial scale as is found \2§\. 

• Diversity and metabolic rate: S"* oc 1 /R so that di- 
versity decreases with increasing metabolic rate: we 
expect ecosystems dominated by endothermic or- 
ganisms with high metabolic rate to have lower di- 
versity than those dominated by ectothermic, low- 
metabolism, organisms (e.g. [1]). 

• Latitudinal gradient rule: For fixed metabolic rate, 
habitat size, and so on, diversity will increase with 
resource. Provided that other factors that link re- 
source uptake to available sunlight do not vary 
our general macroecological relationship encapsu- 
lates the latitudinal gradient rule. Again, this trend 
is present alongside patterns in the other variables 
in ((T3)) as discussed by [l^. For example, p3)) pre- 
dicts that low metabolic rate ecosystems where the 
rate of resource supply is high will be more diverse 
that high metabolic rate ecosystems where the rate 
of resource supply is low; this may suggest a refined 
version of the latitudinal gradient rule. 

• Species Area Relationships: A corollary of length- 
scale dependence of the density is that 5'* oc 
i^~^/G*(^) that is, diversity will vary both 
with habitat size (which is a function of L) and 
the lengthscale of the observation or characteristic 
lengthscale of some process lengthscale (which is a 
function of r*); these are known as Species Area 
Relationships (SAR) (see eg [s^l)- Thus if the in- 
dividuals grow in clumps, say by division, or live 
on a fractal structure (tree, coral, mountain, river) 
or forage in a random walk pattern (jl3|) will con- 
strain the resulting ecosystem SAR. For example 
power law SAR arise if available productive surface 
area or volume orders the availability and uptake 
of resource [31| and that this is in turn ordered by 
the roughness of the terrain which can be modeled 
simply as a fractal [s^, HI] , see also [s^l • Within a 
given ecosystem, diversity and abundance will both 
vary with the r* over which they are observed, as 
well as with each other. If all other variables are 
not strongly varying, their functional dependence 
on can be obtained from the data by the method 
of scaling collapse as has been done by |35| . Again, 
these SAR and the underlying constraint on dis- 
persal and clumping from which they originate will 
also be found to interact with other variables such 
as productivity [s^ and bodysize [s^l- Observa- 
tions over the largest regional or continental scales 



tend to integrate or aggregate over detailed spatial 
dependence and over other variable factors such 
as metabolic rate. These scales exceed that over 
which the terrain varies, and over which processes 
occur that yield spatial clumping, over such large 
scales the effective G — > 1 and the landscape is 
essentially 'flat' so that D — 2. Hence on these 
largest scales the spatial dependence in van- 
ishes, we have that P ~ 5*, and positive relation- 
ships between productivity and diversity emerge as 
has been found 1281. 



Abundance, the 'more individuals' hypothesis: The 
number of individuals in a species n^L^ oc 
aPL? /R so that abundance (density of individu- 
als in each species) increases with productivity in- 
tegrated over the habitat and decreases with the 
typical metabolic rate (e.g. Q)- As we increase the 
total ecosystem energy take up rate, from equation 
(jl3p both the diversity and abundance can increase 
(for fixed metabolic rate). We will find the above 
patterns when one effect does not dominate, for ex- 
ample, Wright's rule will not be seen if the increase 
is entirely in abundance, and not in diversity. The 
relationship between the number of individuals and 
the number of species has long intrigued ecologists, 
and whilst a positive relationship is sometimes as- 
sumed (the 'more individuals' hypothesis for the in- 
crease in richness with overall abundance: see (27| ) 
the mcchanism(s) involved are far from clear. 



• Abundance, diversity and body size: oc l/(n,S',) 
so that the average abundance ^38r,:42] and diversity 
[3] are found to have statistical trends that decrease 
with increasing average body size. 



• // body size trends correspond to trends in trophic 
level then since i?, oc PL^/R the number of 
trophic levels that the ecosystem can support 
increases with total net productivity summed 
over the habitat and decreases with the typical 
metabolic rate [4|. Alternatively the number of 
trophic levels for example may be determined by 
the level of complexity 5* which increases with to- 
tal net productivity summed over the habitat and 
decreases with the typical metabolic rate. 



Our result, equation (jlSp . identifies a functional rela- 
tionship between ecosystem variables. It also immedi- 
ately suggests why many such patterns show consider- 
able variability; for example, plotting abundance versus 
body size does not eliminate other dependencies such as 
on metabolic rate, diversity and the rate of supply of re- 
source. We now use to identify a method to isolate 
these patterns. 
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IV. NORMALIZING TO TEST FOR PATTERNS 
WITHIN AND ACROSS ECOSYSTEMS 

As a starting point say we have a set of observations 
based on individuals consistently organized into cate- 
gories. For each category we observe on lengthscale r* 
the density n,, diversity S"* and body size -B*, so that 
for many such categories we have a set of observations 
of n*, [ni,n2, ■■rij...], of S^,, [Si, S2, -Sj..], and of B^,, 
[Bi, B2, .-Bj..] where each and refer to size i?* 
of the p*'' category. We first consider the case where 
these are all drawn from the same ecosystem and ob- 
served on the same lengthscale so that we have the same 
a, P, L, R, G and r*. There will be a single function ^'(p*) 
which corresponds to this set of observations at different 

. We can write 

^{n,S,B,) - ^^-^^^{n'S'B') 

so that relative to any particular category n'S'B' we can 
obtain ^'(n*S'*_B*) to within a constant. Alternatively 
one can treat as a parametric coordinate, and obtain 
the trend with any one of , 5* , . 

This procedure is illustrated in Figure 1 where we have 
modelled synthetic data. Our synthetic data are gen- 
erated such that there are trends in abundance, diver- 
sity and body size, and also random scatter in all vari- 
ables. Importantly this random scatter is generated to 
be constrained such that all the sampled categories of 
individuals share the same function Vt; the plots show 
that this constraint only emerges in a plot of l/S:,n^B^, 
versus body size i?*. This also offers a method to com- 
pare different categorization schemes which for the same 
ecosystem could yield different 4'(p); one could then in 
principle normalize for (ie compensate for) any 'bias' in- 
troduced by a particular choice of categorization. 

Practically, data is obtained across ecosystems as well 
as within them so that now a, P and L may vary. A 
particular example of this is testing for SAR. To test for 
a SAR one should plot the normalized diversity S: 

^ _ S,n,RB,^ip,) L^-° 

a.P ~ G{r.,L) ^^^> 

versus L with ^(p*) found as above, or a comparison 
made such that ^(p*) is not strongly varying ('similar' 
ecosystems). 

We illustrate this process in Figure 2 where we have 
modeled synthetic data for a species-area comparison. 
We have generated synthetic data in the same manner 
as in Figure 1, such that there are trends in abundance, 
species richness, and body size, and also random scatter 
in all variables, constrained such that all the sampled 
categories of individuals share the same function ^P. In 
addition each group of data is from a different habitat 
size and has a dependence on area predicted by our result 
with power law dependence of G on L. This clear 
pattern in normalized diversity has a corresponding SAR 



pattern of diversity versus area which has considerable 
scatter 



V. CONCLUSIONS 

We have used a 'bottom up' approach to fix the mini- 
mum set of variables needed to specify a generic idealized 
ecosystem. Physical constraints of energy flow and uti- 
lization over the ecosystem, and dimensional constraints, 
then relate the level of complexity that the ecosystem 
can support to intrinsic variables such as density, diver- 
sity and characteristic lengthscales for foraging or dis- 
persal, and extrinsic variables such as habitat size and 
the rate of supply of resource. These constraints hold 
regardless of the details of how a given ecosystem func- 
tions and require only three key assumptions, (i) that the 
ecosystem is in a dynamically balanced steady state, (ii) 
that some (non- unique) method can differentiate differ- 
ent functional groups of individuals within the ecosystem 
and (iii) that these functional groups are connected to 
each other by the flow of resource. We thus find the con- 
straint on the relationships that can exist between these 
(dimensional) ecosystem variables which is reflected in 
observed macroecological patterns. Our result may ex- 
plain why these general, approximate statistical trends 
appear to be so ubiquitous in nature: we obtain these 
patterns without recourse to any detailed information 
about the structure or dynamics of ecosystems or indeed 
how the data is collected. 

Our result fixes the relationship between ecosystem 
variables. This leads to the dimensionless, or normalised 
variables that need to be constructed to isolate the trend 
in one ecosystem variable from another; this provides a 
new method for isolating macroecological patterns. Com- 
parisons could thus be made between datasets by con- 
trolling for (normalizing against) characteristic metabolic 
rate, abundance and diversity in order to isolate the sta- 
tistical pattern with respect to one of these variables. In 
particular this method isolates a function that expresses 
how complex the ecosystem is and it would be intriguing 
to order the data in this way to determine the level of 
complexity of ecosystems that are found in nature, and 
to what conditions they correspond. An example would 
be comparisons across extinct ecosystems, or between ex- 
tinct and contemporary ecosystems, provided a compa- 
rable sample group could be identified. The fact that 
Wright's rule, species area rules and latitudinal gradient 
rules emerge often, but not always, from the observa- 
tional data gathered across ecosystems may reflect vary- 
ing levels of complexity in these ecosystems, or the effect 
of different schemes for categorizing individuals within 
ecosystems and our results provide a method to control 
for this. 

Departures from these statistical patterns where 
ecosystems are similar, and consistently sampled, then 
imply that the system is in a state of rapid change, i.e., 
abundance or diversity explosion or collapse. Any ecosys- 
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tern which is dynamically balanced in the sense discussed 
above will fall within these macroecological patterns, it 
does not need to be a climax or maximum energy utiliza- 
tion system but simply needs to balance the rate of en- 
ergy uptake with that of usage integrated over the ecosys- 
tem. 

Finally, we have identified a dimensionless control 
parameter for ecosystem density and diversity, namely 
(productivity) x (habitat size) /(typical metabolic rate) 
which quite generally emerges from our dimensional anal- 
ysis. This we suggest will be a control parameter in dy- 
namical models for ecosystems based on energy flow and 
conservation and will order the emergent behaviour of 
these models. We relate this control parameter to the 
level of complexity that a dynamically balanced ecosys- 
tem can support (its order parameter). This control pa- 
rameter is the ratio of energy input rate to the ecosystem 



to the metabolic rate of the smallest possible unit of life, 
a single cell. If it is reasonable to identify a smallest pos- 
sible unit of life then this parameterizes the threshold at 
which life, defined in this manner, can occur (c.f. Love- 
lock's conjecture about dimensionless control parameters 
for prebiotic planets [43l|). 
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FIG. 2: Log-log plot of synthetic data generated for cate- 
gories of individuals from habitats of different areas. The 
data is constructed with power law dependence of diversity 
on lengthscale and trends in abundance, diversity, body size 
and area with random scatter. This pattern is only revealed 
in a plot of dimensionless diversity versus area. 



FIG. 1: Log-log plots of synthetic data generated for an 
ecosystem where all the observed individuals are constrained 
to share the same ecosystem complexity function ^. The 
data is constructed to show trends in abundance, diversity 
and body size with random scatter. The constraint can only 
be discerned by plotting l/SjUjBj versus body size Bj. 



